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ABSTRACT 






Aims. The combined point spread function (PSF) of the Broadband Filter Imager (BFI) and the Solar Optical Telescope (SOT) onboard 
the Hinode spacecraft is investigated. 

Methods. Observations of the Mercury transit from November 2006 and the solar eclipse(s) from 2007 are used to determine the PSFs 

C^^ I of SOT for the blue, green, and red continuum channels of the BFI. For each channel large grids of theoretical point spread functions 

are calculated by convolution of the ideal diffraction-limited PSF and Voigt profiles. These PSFs are applied to artificial images of an 

eclipse and a Mercury transit. The comparison of the resulting artificial intensity profiles across the terminator and the corresponding 

observed profiles yields a quality measure for each case. The optimum PSF for each observed image is indicated by the best fit. 

Results. The observed images of the Mercury transit and the eclipses exhibit a clear proportional relation between the residual 

Q . intensity and the overall light level in the telescope. In addition there is a anisotropic stray-light contribution. These two factors make 

Vh ' it very difficult to pin down a single unique PSF that can account for all observational conditions. Nevertheless the range of possible 

Ti ' PSF models can be limited by using additional constraints like the pre-flight measurements of the Strehl ratio. 

j^ I Conclusions. BFI/SOT operate close to the diffraction limit and have only a rather small stray-light contribution. The FWHM of the 

PSF is broadened by only ~ 1 % with respect to the diffraction-limited case, while the overall Strehl ratio is ~ 0.8. In view of the 

large variations - best seen in the residual intensities of eclipse images - and the dependence on the overall light level and position in 

CO ' the FOV, a range of PSFs should be considered instead of a single PSF per wavelength. The individual PSFs of that range allow then 

^ I the determination of error margins for the quantity under investigation. Nevertheless the stray-light contributions are here found to be 

>0 , best matched with Voigt functions with the parameters cr = 0'.'008 and y = 0'.'004, 0'.'005, and 0'.'006 for the blue, green, and red 

CO ■ continuum channels, respectively. 

lO ■ 

^^ ■ Key words. Sun: atmosphere; Instrumentation: high angular resolution; Space vehicles: instruments 

o 

QQ 1 . Introduction and scattering. Seeing has mostly an impact on the central parts 
f~^ . ^ ... . _ , . . , of the PSF, blurring the image, but also contributes to the wings, 
^ , The quantitative interpretation of an observation requires a good wher eas scat tering becomes significant for the far wings of the 
. ^ ■ knowledge of the properties of the utilised instrument as it can p^p ,JMEtiaQ983). It depends on the elevation of the line of 
K^ Significantly aff^ect empirically derived properties. A realistic ^^ j^^ ^^^ ^^^^^^^ conditions (e.g. dust, aerosols or water vapour 
iVj model or the instrumental influence is or particular importance • j. atmnsnhprp't 
j^ ■ for the comparison of certain aspects of observations and numer- 
■ - - ' ical models, e.g . for die intensity contrast of solar granulation Space-borne instruments have the advantage of not being 
dNordlundll 1 984) limited by the Earth's atmosphere but still instrumental stray- 
Each optical instrument causes the radiation emitted by a light has to be considered. Often the stray-light contribution is 
point source (or an object) to be spread on the image plane anisotropic and varies over the field of view (FOV), making it 
of the detector. The image of an extended light source like the hard to derive a detailed model for the PSF This effect is akeady 
Sun is degraded this way and small spatial scales remain unre- visible in solar limb observations, where the intensity does not 
solved. Mathematically this phenomenon is described by means decrease to zero but to a finite valu^J The instrumental devia- 
of a point spread function (PSF) or its Fourier transform, the tions from the diff^raction-limit are mosdy caused by imperfec- 
Optical Transfer Function (OTF). Even an otherwise perfect op- tions and contamination of the optical components in the tele- 
tical telescope induces image degradation due to diff^raction at scope, e.g. dust on optical surfaces or impurity of lens material 
the finite aperture stop. The resulting diffraction-limit eff^ectively but also microscopic scratches or micro-roughness. Also subde 
defines the best theoreticaUy achievable angular resolution of a eff^ects like jitter of the spacecraft can contribute to the PSF The 
telescope resulting (anisotropic) stray-light increases the relative intensity 
In reality there is additional degradation due to (i) the Earth's contribution from the PSF wings and with it the relative height of 
atmosphere (for ground-based instruments) and (ii) instrumental the side-lobes, which nevertheless stay in position. With a higher 
eff-ects. The terrestrial atmosphere introduces image degradation contribution from the wings the relative height of the central 
due to turbulence along the line of sight - the so-called seeing - 
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peak and with it the corresponding Strehl raticQ decrease, while 
the central peak is broadened. 

Constructing a realistic PSF is unfortunately a non-trivial 
task as often important contributions are hardly known. This is 
in particular true for stray-light produced inside an optical in- 
strument (Mattig 1983). The best way to determine the proper- 
ties of a telescope is of course the direct measurement in the 
laboratory but for space-borne instruments it is usually diffi- 
cult to produce the same conditions as for the later operation. 
Another possibility of measuring the PSF is the observation of 
a partial solar ecli p se. A n illustrative example is the work by 
iDeubner & Mattig* (Il975h who used observations of a partial 
eclipse in 1973. They fitted the observed intensity profile across 
the eclipse t erminator with a combination of two Gaussians (cf. 
lLevvll97lh . lNordlundl (119841 ) argues that the wings of the spread 
function give rise to uncertainties, which make it problematic 
to retrieve a definite analytical fit of the measured profile. He 
showed that instead of two Gaussian also a combination of two 
Lorentzians produces a fit of comparable quality. The related un- 
certainties can have a fundamental impact on the derived inten- 
sity contrast. Nordlund proposes a correcti on factor of 1 .5 with 
respect to the results of iDeubner & Mattia Beside these uncer- 
tainties, a fundamental problem remains as the fit of an observed 
terminator profile usually assumes a symmetric PSF/OTF but 
hardly can account for asymmetric contributions. As we shall 
see in the case of SOT a real PSF can have significant anisotropic 
components. 

In many cases, however, suitable observations of an 
eclipse or a transit of Mercu ry or Venus (see, e.g., 
iDeForest & Wills-Davevl Isubmittedl) are not available. Instead 
one has to estimate PSF properties from the spatial power spec- 
tral density and other ad hoc arguments. Another possibility 
is to estimate the PSF by comparing the intensity distribution 
(histogra ms) of a parti cular obse rvation with numerical models 
(see., e.g-.ISchiissler et al. 2003; Leenaarts & Wedemever-Bohml 
120051: iLangangen et al.1,2007,) . Often the PSF is approximated 
with a linear combination of an Airy and a Lorentz func- 
tion, characterised by the two parameters a, the height of the 
Lorentzian relative to the height of the Airy core, and /3, the 
FWHM of the Lorentzian wings. In most cases, however, there is 
an extended range of possible combinations of a and/3. More se- 
vere is the (implicit) assumption that the numerical model, which 
is used as reference for the fit between observed and synthetic 
image, is perfectly realistic. 

H ere a set of PSFs for the Solar Optical Telescope 
(SOT , iTsuneta et al."2008'; 'I chimoto et al]|2008t ISuematsu etal] 
l2008t ISh imizu et al. 2007) onboard the Hinode spacecraft 
dKosugi et al.ll2007) is derived from observations of a Mercury 
transit and two eclipses for different wavelengths. The ob- 
servations are described and analysed in Sect. |2] Theoretical 
two-dimensional PSFs are constructed by convolution of ideal 
diffraction-limited PSFs and Voigt functions and applied to an 
artificial eclipse and a Mercury transit in Sect. [3] By comparison 
of the resulting intensity profiles and the observational results, 
the best PSF is determined for each case. Discussion and conclu- 
sions follow in Sects. |5]and|6l respectively. SOT observations of 
quiet Sun granulation will be systematically analysed and com- 
pared to synthetic intensity images from numerical simulations 
in forthcoming papers. 



^ The Strehl (intensity) ratio is defined as the ratio between the cen- 
tral peak values of a PSF and the one of an ideal, i.e. aberration-free, 
difFraction-limited PSF. It provides a measure for the deviation from 
the diffraction limit. 



2. Observations 

2.1. Filtergrams 

The observations analysed in this study were obtained with the 
Broadband Filter Imager (BFI) of the Solar Optical Telescope 
(SOT) onboard the Hinode satellite (see Sect.[T]for references). 
Primarily the blue, green and red continuum channels of the 
BFI are considered. The central wavelengths of these filters are 
450.45 nm, 555.05 nm, and 668.40 nm, respectively. All three fil- 
ters have a band width of 0.4 nm. In addition a few images for 
the other three BFI filters are used. These filters are CN, Ca II H, 
and G-band with central wavelengths of 388.35 nm, 396.85 nm, 
and 430.50 nm, respectively. The correponding band widths are 
0.7 nm, 0.3 nm, and 0.8 nm. 

Observations of the Mercury transit from November 8th, 
2006, and the solar eclipse of March 19, 2007, are analysed (see 
Fig. [T] for examples). The eclipse was partial on Earth but total 
for Hinode. It is referred to as total eclipse hereafter. There also 
occured a partial eclipse on February 17, 2007, which was ob- 
served by Hinode but which was not visible from the ground. 
It was observed with the G-band channel only but still pro- 
vides useful information for this study. Additional observations 
of quiet Sun granulation for various positions on the disc and at 
the solar limb are used. 

The position of the FOV and individual subregions on the 
solar disc is here reduced to the heliocentric position fi = cos 6. 
No distinction is made between the different parts of the limb 
(N/W/S/E) as there is no noticeable difference in the limb pro- 
files for the quiet Sun (see, e.g. Langh ans & Schmidt 200 2). The 
pixel scale of the detector is '.'05448 and thus smaller than the 
spatial resolution (and even the diffraction limit) of the optical 
system (0'.'2 - 0'.'3, see, e.g., Kosugi et al. 2007, and Sect.[T]). As 
the same is still true for images with 2x2 pixel-binning, no dif- 
ference in terms of contrast etc. is to be expected. Both binned 
and unbinned images are thus used next to each other. 

The data reduction (incl. dark subtraction and flat-fielding) 
was caiTied out with Solar Soft (Freeland et al. 2000). Dark 
frames and flat-fields were kindly provided by T. Berger. In a 
few cases remaining hot pixels, which often occur in compact 
groups of adjacent pixels, had to be removed. 

Residual intensities. The intensity profiles and residual inten- 
sities in the dark areas are determined for the filtergrams from 
the Mercury transit, the two eclipses, and from limb observa- 
tions. The residual intensity includes contributions from the PSF 
wings due to the optical system but also contributions from scat- 
tered light. All results are summarised in Table [1] for different 
distances from the terminator (for the eclipses) and the solar 
limb, respectively. For the Mercury transit, the residual inten- 
sity at its disc-centre is given (r ^ 5"). The table contains the 
values for /res, which are the averages over all profiles, the corre- 
sponding standard variation, and (in brackets) the extrema. A^ is 
the number of considered images. The intensity profiles are aver- 
aged parallel to the terminator/limb, resulting in a radial average 
for the Mercury transit. 

Figure [T] illustrates the determination of the residual light 
level for the case of the green continuum. For the partial eclipse, 
however, the G-band is chosen because no other observations 
are available. The right column displays the intensity profiles 
derived from the images on the left (black solid lines). In ad- 
dition also the corresponding profiles from all other images for 
the same wavelength are plotted as grey lines in each panel. All 
profiles are averaged over pixels with the same distance to the lu- 
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Table 1. Residual intensities /res in the observations at different distances r from the terminator or the solar limb, respectively. See 
text for details. 



band 
A [nm] 



CN 
388.3 



CallH 
396.9 



Gband 
430.5 



be 
450.5 



555.0 



re 
668.4 









total eelipse 






A' 

r = 5" 
r = 10" 
r = 20" 


4 
2.0 ±0.2 [1.7, 2.3] 
1.5 ±0.2 [1.2, 1.7] 
0.8 ±0.6 [0.0, 1.3] 


3 
2.3 ± 0.6 [2.0, 3.0] 
1.9 ±0.5 [1.6,2.4] 
1.5 ±0.4 [1.3,2.0] 


3 
3.4 ±0.9 [2.8,4.4] 
2.9 ±0.7 [2.4,3.7] 
2.6 ±0.6 [2.0,3.2] 


3 
2.9 ±0.4 [2.4, 3.2] 
2.3 ±0.3 [2.0, 2.5] 
1.9 ±0.2 [1.7, 2.1] 


8 

4.7 ±1.5 [1.7,5.9] 
4.2 ±1.5 [1.1,5.3] 

3.8 ±1.4 [0.8,4.8] 


3 

5.3 ±1.5 [3.6,6.4] 
4.7 ± 1.5 [3.0, 5.7] 

4.4 ±1.4 [2.8,5.2] 



partial eelipse 



A' 






9 








r = 5" 






2.0 ±0.8 [1.1, 3.2] 








r = 10" 






1.6 ±0.7 [0.9, 2.7] 








r = 20" 






1.4 ±0.6 [0.7, 2.2] 









Mereury transit at solar disc-centre 



A' 
r = 5" 



5.9 ± 0.9 [4.9, 6.7] 



6.2 ±0.3 [5.7, 6.5] 



8.1 ±0.9 [7.1,9.2] 



7.0±0.7 [6.2,8.1] 



10.3 ± 0.2 [10.2, 10.6] 11.3 ± 0.5 [10.7, 12.0] 



Mercury transit at solar limb 



A' 

r = 5" 



39 
8.0 ±2.2 [6.3,13.: 



solar limb 



A' 








66 


18 


20 


r = 5" 








5.3 ±0.7 [4.3, 6.8] 


6.9 ±1.1 [5.7,10.5] 


7.9 ±1.3 [5.5,9.9] 


r = 10" 








4.7 ±0.6 [3.7, 6.1] 


6A±1.2 [5.1,9.9] 


7.2 ±1.2 [5.0,9.1] 


r = 20" 








4.0 ± 0.7 [2.9, 5.7] 


6.0 ±1.3 [3.6,9.0] 


6.7 ±1.1 [5.0,8.3] 



nar/solar/Mercury limb, respectively. There is a general trend of 
a higher residual intensity with increasing wavelength, although 
the mean G-band value is slightly larger than for the neighbour- 
ing bands (see Table [TJ. The variation of /,e,s,,i for images of the 
same wavelength is significant, clearly demonstrating that it is 
impossible to describe the stray-light contribution with a single 
number (for each wavelength). The situation is much better for 
the Mercury transit at solar disc-centre as the individual images 
were taken under very similar conditions, resulting in less vari- 
ation of /les- In the following the individual sets of images are 
characterised. 



Total eclipse. The mean residual intensities in the images for 
the total solar eclipse at a distance of 5 "from the terminator 
amount to (2.9 + 0.4) % of the mean intensity of the unocculted 
solar disc for the blue continuum but to (5.3 ± 1.5) % for the red 
continuum, respectively. At a distance of 20 "the values reduce 
to (1.9 + 0.2)% and (4.4 ± 1.4)%. Values between 0.8% and 
5.2 % are found for individual (averaged) profiles. 

Partial eclipse. For the partial eclipse only G-band images are 
available. The intensity profiles essentially form two groups, one 
with a residual intensity at a distance of 20 "of 0.7 % to 1.0 % 
of the solar disc intensity and the other with 1.8 to 2.2 %. The 
higher values result from images where the dark occulted area is 
less than 42 % of the image, resulting in a higher mean count rate 
(see below). The low values come from images where the lunar 
disc covers more than half of the image and the mean count rate 
is low. 



Mercury transit at solar disc-centre. Most images of the 
Mercury transit on the solar disc are taken at a heliocentric po- 
sition of the Mercury disc of yU a; 0.9. The residual intensity 
at Mercury disc-centre is determined from the radially averaged 
profile across the Mercury limbwith respect to the mean of sur- 
rounding solar disc area of comparable ju. Close to solar disc- 



centre (ji a: 0.9) the following values are derived for the three 
continuum bands: (7.0 ± 0.7)% (blue), (10.3 ± 0.2)% (green), 
and (11.3 + 0.5) % (red). For the other three broadband channels 
the following values are found: CN (388 nm): (5.9 ± 0.9) %, Ca 
(397 nm): (6.2 ± 0.3)%, and Gb (431 nm): (8.1 ± 0.9)%. The 
Mercury transit at solar disc-centre obviously gives residual in- 
tensities that are much higher than estimated from the eclipses 
and from the solar limb observations. This can already be seen 
from the intensity profiles in Fig.[Tl which decrease to much 
lower residual intensity values off the solar limb than compared 
to the centre of the Mercury disc. This is partly due the influ- 
ence of the wings of the telescope's PSF, which contribute more 
to the residual intensity for the Mercury disc than for an eclipse 
(see Sect.m. 



Mercury transit at the solar limb. In addition there are green 
continuum images at the solar limb covering the beginning of the 
Mercury transit. The value of /,es,,i referres to the mean intensity 
of the solar disc close to the Mercury disc in order to ensure a 
comparable /i. Next to the radial profiles also profiles along the 
X- and y-axis through the Mercury disc-centre are shown in the 
left part of Fig.[l]. All include also the intensity decrease over 
the solar limb. At Mercury disc-centre the radial averaged pro- 
files exhibit residual intensities of (8.0 ± 2.2)% with extreme 
values of 6.3 % and 13.8 %. Only the inner parts, which are of 
interest here, of the radial profiles are reliable as further away 
from the disc-centre one averages the non-radial solar disc. The 
averages along the x- and y-axes show very similar residual in- 
tensities (see Fig.[T|). Nevertheless the derived values for /res,^ 
are smaller than for the Mercury transit at disc-centre because 
now a larger area of the image is dark, mapping the space be- 
yond the solar limb. The large variation in /les,^ can be explained 
by (i) the varying fraction of the dark area and thus the overall 
light level in the telescope and (ii) the more problematic deter- 
mination of a reference intensity as the Mercury disc at the solar 
limb spans a larger range in //. 
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Solar limb. The pure limb observations (i.e. without Mercury, 
see Fig. [T^-b) imply values for the residual intensity from 
(5.3 ± 0.7)% for blue continuum to (7.9 + 1.3)% for the red 
continuum at a distance of 5 "from the solar limb. For 20 "the 
values decrease to (4.0 ± 0.7) % and (6.7 + 1.1) %, respectively. 
The residual intensity refers to the mean intensity of the solar 
disc close to the limb (average running from -1 "to -2"). The 
reference, however, is problematic as the limb images often ex- 
tend to // ~ 0.4 and in some cases even to yu > 0.5. 

2.2. Dependence on overall light level 

The large variation in the eclipse intensity profiles for the same 
wavelength is caused by the different illumination of the detec- 
tor during the exposure. This "overall light level" is here quan- 
tified as the "mean count rate" C. It is defined as the total count 
number of all pixels, i.e. the full FOV, of an image divided by 
the FOV area (in arcsec^) and the exposure time. For a given 
wavelength, an unocculted image taken at solar disc-centre has 
a higher value of C than an image taken closer to the solar limb. 
C is reciprocally proportional to the fraction of the occulted and 
thus dark image area. For the images considered in this section, 
C is highest for the Mercury transit at disc-centre. In all other 
cases a significant part of the detector is dark as it images the 
lunar disc or the background beyond the solar limb. The same is 
true for the Mercury transit images at the solar limb. 

The residual intensities are measured for all images at a dis- 
tance of 5" from the terminator/solar limb. The results for the 
different kinds of observations are represented by different sym- 
bols in Fig. |2] For all six wavelength channels I^es.A increases 
linearly with C. The eclipse and Mercury transit data was used 
to derive the regression lines shown in Fig. |2] The solid lines 
show the result of a linear regression, which reproduce the rela- 
tion quite well. The solar limb profiles are excluded here because 
of the uncertain reference intensity. The solar limb data is in line 
with the eclipse and Mercury transit when plotting the corre- 
sponding count rate instead of /res. Obviously a higher reference 
intensity is necessary to make the limb data comply with the re- 
gression lines for /res. The dependence on C indicates that the 
residual intensity contains contributions from anisotropic stray- 
light, which, e.g. may originate from internal reflections in the 
instrument and thus clearly depend on the overall light level. 

The heliocentric positions fi are plotted as function of 
the mean count rate for the blue, green, and red continuum 
in the panels e, h, and k of Fig. |2] respectively. The pan- 
els f, i, and 1 show the residual intensity as function of fi. 
The Mercury transit images for the green continuum have 
residual intensities of 2.6% to 4.8% close to the solar limb 
(position of the Mercury disc at /i < 0.2), whereas ~ 10.0 % are 
derived close to disc-centre. It clearly demonstrates the increase 
of the relative residual intensity with fi, related to the centre-to- 
limb variation of the mean intensity of solar granulation and the 
corresponding change of the mean count rate. 

The residual intensities in images of "regular" unocculted 
granulation are now estimated from the corresponding values of 
C and the (wavelength-dependent) relation with /res,,;, which was 
found from the Mercury transit and the eclipse data. The results 
for the blue, green, and red continuum are represented as dots in 
Fig.|2l Finally a linear regression of the form I^e^j = m ■ /u + c 
is performed for each wavelength (solid lines in Fig.|2l). The de- 
rived parameters c (constant) and m (slope) are listed together 
with the maximum value of /res,^ at // = 1.0 for all three consid- 
ered continuum wavelengths in Table |2] At solar disc-centre the 
following residual intensities Ii-g^.A can be expected in images of 




3.0 -2.5 -2.0 
X [10^ arcsec] 



Fig. 3. Stray-light in a SOT green continuum filtergram of the 
total eclipse from March 19th, 2007 (0.78 < fi < 0.92): Image 
(left) and profiles along the x-axis (right). In order to make the 
stray-light visible, the data range is limited to a value of 120, 
which is about 10 % of the mean value of the bright solar disc 
(lower left corner of the image). The vertical positions of the pro- 
files are indicated by tick marks on the right side of the image. 
Profile a lies in a "shadow", whereas profile b crosses a "streak" 
at a position of x = -300 ". 



"regular" granulation: 6.4 %, 9.5 %, and 10.9 %, for blue, green, 
and red continuum, respectively. The estimate is based on resid- 
ual intensities at 5 "distance from the terminator of the Mercury 
transit and the eclipses and should thus be interpreted as upper 
limit only. 



2.3. Anisotropy 

A direct proof of stray-light is presented in Fig. |3] It shows a 
frame taken during the total eclipse but with the displayed inten- 
sity range limited to only 10% of the mean value of the bright 
solar disc. There are reflections (e.g. from baffles) and a shadow 
at the upper edge of the detector More "artefacts" that are not 
removed during the data reduction process can be found. For 
instance, there is a streak-like feature close to the terminator. 
Obviously the stray-light is highly anisotropic. All those details 
are relatively subtle and might be negligible for most applica- 
tions for which exact photometry is not needed. Nevertheless the 
residual intensity can reach up to 4 % of the solar disc intensity 
in this example. Similar elTects are expected for the other chan- 
nels, too. For instance, the red continuum channel of the BFI 
shows a prominent spot. Variations in the stray-light contribu- 
tion due to asymmetric scattering inside the o ptical system have 
been reported for other telescopes before (e.g. lMattial 19831) . 



3. Theoretical PSFs for SOT 

3. 1 . Ideal diffraction-limited PSFs 

The diffraction-limited PSF of a cuxular homogenous illumi- 
nated aperture of diameter d is described by the Airy function. 



Table 2. Estimate of the residual intensity I^es.A for images of 
"regular" granulation in percent points. See text for details. 

"Wavelength m [%] c [%] (/,.es.i).,ax [%]~ 



450.45 
555.00 
668.40 



3.18 
3.77 
4.51 



3.22 
5.73 
6.39 



6.40 
9.50 
10.90 
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Fig. 4. Theoretical point spread functions for a wavelength of 
555 nm: Original diffraction-limited PSF (top row), after con- 
volution with a Voigt function with y - 0'.'004 and cr - 0'.'012 
(middle), and after linear combination with a Lorentzian with the 
parameters a - 0.25 and /3 = 0.05 (bottom). The left columns 
shows the inner part of the PSFs on a logarithmic grey-scale, 
which covers six orders of magnitude, whereas the correspond- 
ing profiles along the x-axis are plotted in the middle column. 
The rightmost panels display the the PSF core on a linear scale 
(solid) together with the original PSF (dotted). The dashed line 
in the lower right panel is the used Lorentzian profile. 



For a wavelength A the FWHM of the central Airy disc is equal 
to the resolution of the telescopqj 



smares - - 



(1) 



For a more complicated aperture geometry the PSF can be 
calculated as the square of the wave number Fourier transform 
of the pupil shape. It has to be normalised so that the sum is 
equal to one. For SOT the 50 cm diameter aperture, the 17.2 cm 
diameter central obscuration, and the three spiders (120° sepa- 
ration, 4 cm wide) are taken into account. Each PSF is defined 
as an array with 5120^ pixels and a spatial increment of half a 
SOT detector pixel. The extent of the PSFs is thus larger than the 
minor axis of the FOV of BFI images and includes the far wings. 
The central lobes have the following FWHMs (cf. Eq. ©): O'.'IS 
at 450.45 nm (blue), 0'.'22 at 555.05 nm (gree n), and 0'.'27 at 
668.40nm (red), respectively. [Suematsa (l2007h kindly provided 
an ideal G-band PSF for the main mirror of SOT, which has a 
FWHM of O'.'IS. As ideal PSFs can be scaled reciprocally with 
wavelength, the G-band PSF can be compared to the PSFs cal- 
culated here. They agree within the limits of numerical uncer- 
tainties. 



^ It should not be confused with the Rayleigh criterion, which is de- 
fined with an additional factor of 1.22 for the position of the first min- 
imum surrounding the central Airy disc (approximately the first mini- 
mum of the Bessel function of the first kind). 



3.2. Non-ideal PSFs 

Here Voigt functions are used as approximation for the non-ideal 
PSF contributions (see Sect.!!). The final PSFs are then derived 
by a convolution of the ideal PSFs Pn and a Voigt function 'Vy 



%,,x(r) = (2 






+O0 g-r'V(2o-^) 



dr\ 



which is the convolution of a Gaussian 

g-'-/(2o--) 

V27ro- 
and a Lorentz function 



-Cy(r)^ 



n(r^+'-) 



(2) 



(3) 



(4) 



Here r - -\Jx^ + y^ is the spatial ordinate and y is the FWHM 
of the Lorentzian. The FWHM of the Gaussian is defined as 
2 cr(2 In 2)'^^. With u - r/( ^a) and the damping parameter 



r 

a — — ■ 



1 



2 ^|2o- 
the Voigt function can also be written in the form 



'V{a, u) = 



I a p 

V27rcr ^ J-a 



a^ + (u — u') 



■ du — 



H(a, u) 



2na- 



(5) 



(6) 



The far wings, i.e. the PSF at large angles, are normally domi- 
nated by the Voigt function. The function is numerically reduced 
to a delta peak if the FWHM of the Voigt function is signifi- 
cantly smaller than the numerical axis increment. In that case 
the convoluted PSF is identical to the original diffraction-limited 
one. Here PSFs with an additional Voigt function 'Vy.a- are com- 
puted for a parameter grid of y and cr. The grid includes PSFs for 
which Vdx is convolved with a pure Lorentzian Xy for different 
values of y or with a pure Gaussian Qcr for different values of cr. 
In some studies a linear combination of the ideal PSF and a 
Lorentzian of the type 



^L (-«, y) = ( 1 - a) ^di (x, y) + aj:p (r(x, y)) 



(7) 



is used instead of a convolution. The linear approach is consid- 
ered here only for comparison as it is strictly speaking physi- 
cally not correct. Instead the convolution method is preferrable 
(see Sect.|5]l. To avoid confusion with the convolution approach 
the FHWM of the Lorentzian is here named yS instead of y. A 
parameter grid for [a,li\ is computed. 



4. Fitting observations with synthetic profiles 

Here the different theoretical PSFs for SOT (see Sect. |3j are ap- 
plied to artificial idealised two-dimensional images of a Mercury 
transit and an eclipse. The image dimensions conform to half 
of the BFI detector, i.e. 2048 x 2048 pixels with a scale of 
'.'05448 px arcsec-i , resulting in a field of 11 1 '.'6 x 11 1 '.'6. The 
resulting synthetic degraded intensity profiles are then compared 
to empirical profiles derived from SOT images (see Fig.|5]for an 
example). The optimum PSF is the one that produces the small- 
est difference between observed and synthetic intensity profile. 
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Fig. 5. Artificial observation of a Mercury transit at a wavelength 
of 555 nm: a) innermost region of the original intensity mask, 
b) degraded image after convolution with a non-ideal PSF that 
was calculated by a diffraction-limited PSF and a Voigt profile 
with y - 0'.'012 and cr - 0'.'098, c) radially averaged intensity 
profiles for the original image (dotted) and the degraded image 
(solid), and d) close-up of panel c. The dot-dashed line is an 
empirical intensity profile derived from a SOT image. 



4.1. Ideal diffraction-limited PSFs 

4.1.1. Artificial Mercury transit 

The ideal theoretical PSFs for each wavelength are applied to a 
simple binary mask of a bright solar background with a black 
Mercury disc (cf. Fig.|5]l. Convolution with the PSF transforms 
the sharp edge of the Mercury limb in the initial mask into a 
smooth slope. The ideal SOT PSF for a wavelength of 450.1 nm 
(blue continuum) produces a residual intensity /j-es of 1.7% 
(with respect to the bright solar background) at the centre of 



Table 3. Residual intensities /,es in the artificial diffraction- 
limited^ and the SOT observations'". 



the Mercury disc in the degraded image. A simpler Airy PSF 
results in only 0.7 % residual intensity, implying a substantial 
contribution due to the spiders and the central obstruction of the 
telescope. Already replacing the SOT PSF with its radial mean, 
i.e., removing the non-radial variation due to the spiders, signifi- 
cantly decreases the residual intensity. Quantitative results obvi- 
ously demand for a detailed model of the telescope properties. 

Although the considered FOV is already relatively large, the 
residual intensity still grows asymptotically as a function of im- 
age (and PSF) size. The far wings of the PSF contribute only 
little to the intensity but nevertheless produce an asymptotic in- 
crease of /les (at disc-centre) by a factor of ~ 1.1 with respect to 
a field size of 1 1 1 .6 " x 1 1 1 .6 ". Approximately the same factor 
is found for A - 555.0 nm and 668.4 nm. The estimate is based 
on an analytical 2D Airy function. Under the assumption that the 
detailed PSFs show the same qualitative behaviour, the factor 1 . 1 
is now multiplied to the corresponding values of /,es, resulting in 
residual intensities of 1.8% (450.1 nm), 2.2% (555.0nm), and 
2.6% (668.4 nm), respectively (see Table |3|l. The proportional 
increase with wavelength is related to the broadening of the PSF 
(see the FWHM). 



4.1.2. Artificial eclipse 

The experiment is repeated for artificial eclipse images with a 
straight terminator Different field sizes and also positions of the 
limb within the images are considered. Again the residual in- 
tensity, which is derived with the ideal PSF, is corrected for the 
asymptotic behaviour for very large bright areas due to the far 
wings. At a distance of 5 ", which corresponds to the Mercury 
disc radius, the following residual intensities (corrected with a 
factor of 1.1 for the far wings) are found to be 0.6%, 0.7%, 
and 0.8 % for blue, green, and red continuum, respectively (see 
Table O. The values are much smaller than for the Mercury ex- 
periment owing to the different geometrical situations and the 
different proportions of bright to occulted area. In contrast to 
the Mercury transit experiment, where the geometry is radial- 
symmetric and the occulted area is relatively small, the eclipse 
has a clear anisotropic geometry with a significant occulted area. 
At a location beyond the eclipse terminator only radiation from 
one side is received, whereas the centre of the Mercury disc is 
irradiated by the bright (solar) disc from all directions. 



/res[%] 



450 nm 



555 nm 



668 nm 



Mercury disc-centre (r = 5 ") 



r = 5" 


observed 


7.0 ±0.7 


10.3 ± 0.2 


11.3 ±0.5 




theoret. 


1.8 


2.2 


2.6 




difference 


5.2 ± 0.7 


8.1 ±0.2 


8.7 ±0.5 



total eclipse, different distances r from lunar limb 



r = 5" 


observed 


2.9 ± 0.4 


4.7 ± 1.5 


5.3 ± 1.5 




theoret. 


0.6 


0.7 


0.8 




difference 


2.3 ± 0.4 


4.0 ± 1.5 


4.5 ± 1.5 


;•= 10" 


observed 


2.3 ± 0.3 


4.2 ± 1.5 


4.7 ± 1.5 




theoret. 


0.3 


0.3 


04 




difference 


2.0 ±0.3 


3.9 ± 1.5 


4.3 ± 1.5 


r = 20" 


observed 


1.9 ±0.2 


3.8 ± 1.4 


4.4 ± 1.4 




theoret. 


0.2 


0.2 


0.3 




difference 


1.7 ±0.2 


3.6 ±1.4 


4.1 ± 1.4 



^ All theoretical values are corrected for the influence of the far PSF 
wings at different distances from the terminator (see Sect. l4.1.Tt . 

'' The difference between both is due to contributions not included in 
the ideal PSFs, e.g., stray-light inside the telescope. 



4.1.3. Comparison with observations 

Comparison of the observed and artificial intensity profiles show 
that when using an ideal PSF (i) the average observed intensity 
profiles over the terminator have a less steep gradient and (ii) the 
residual intensities derived from the observations are generally 
larger (see Table O. This indicates that (i) the true PSFs are 
broader than the theoretical diffraction-limited ones and (ii) that 
there is an additional contribution due to stray-light. The differ- 
ence between theoretical and observed profiles decreases with 
distance from the terminator and can thus not be described with 
just a constant intensity offset. Many profiles, averaged parallel 
to the terminator or just individual profiles, exhibit even a change 
in slope. The slope can get steeper but also reverse again in some 
cases (see Figs. [T]and|3]l. The synthetic degraded intensity pro- 
files would fit the terminator slope of the observed (averaged) 
profiles better if the PSFs were broader A broader PSF corre- 
sponds to a smaller effective aperture (or a longer wavelength, 
see Eq. ([T]i). However, even though a broader PSF would provide 
a better fit of the slope, the (dark) part of the profiles directly be- 
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yond the terminator is often not well reproduced. Many observed 
profiles do not show the theoretical asymptotic decrease towards 
a residual intensity, indicating that a purely diffraction-limited 
PSF is too simple. 

4.2. Non-ideal PSFs incl. convolution with a Voigt function 

Now those PSFs are used that are calculated by convolution of 
the diffraction-Umited PSFs and a Voigt function. Synthetic in- 
tensity profiles are again derived from degraded images of the 
artificial eclipse and the Mercury transit and are then compared 
to the observed intensity profiles. An example for the Mercury 
transit is displayed in Fig|5] 

A constant intensity offset Iq, which is added to the synthetic 
profiles /syn.orgi is introduced as a third parameter: 



/syn (X, y, Cr,Io)^h + {^- h) /syn,oi-g {x, J, Cr) 



(8) 



The discrepancy between observed and synthetic profiles is mea- 
sured for all parameter combinations [y, cr, /o]. The quality of a 
fit is evaluated with the quantity 



E (y, cr, /o) = 



lobsix) - hya{x,y,cr,Io) 

/syn {X, y, cr, /o) 



(9) 



averaged over the terminator region and the first 10 "of the dark 
part in case of an eclipse and over the first 5 "in case of the 
Mercury transit), respectively. This goodness-of-fit measure is in 
principle a sum-of-squares or least squares approach. For every 
[y, cr] the fitting procedure returns an error vector for the dif- 
ferent /q. Although the general picture remains the same when 
choosing a fixed !„, the minimum error is chosen for each [y, cr], 
effectively reducing the measured errors to Sminiy, cr)- The pro- 
cedure is repeated for every observed profile. 

Reference case. First a reference case is constructed in order 
to investigate the fitting routine and the quality of the selected 
error measure. Instead of fitting an observed intensity profile, a 
known synthetic profile for y - O'.'OOS and cr - O'.'Ol is preset. 
The error £„,!„ forms a distinct pattern in the y-cr plane, which is 
shown in Fig.|6] The black track, where the error is smallest, is 
called here the "minimum error trench" (MET). The properties 
of the Voigt function make the MET bend and form a "knee". 
In the reference case depicted in Fig. |6] the "knee" is located at 
y X 0'.'03 and cr » 0'.'016. Below the knee, i.e. for smaller cr, 
the MET follows lines of constant Strehl but also roughly lines of 
equal broadening of the FWHM of the PSF core with respect to 
the pure diffraction-limited PSF. In contrast the part of the MET 
above the knee follows a line of constant damping parameter a 
of the Voigt function (see Eq. (|5]l). The MET in the reference 
case fades again along the upper right part, indicating that the fit 
between the intensity profile degenerates for increasing cr. The 
determined best-fitting parameter pair (marked with a circle in 
the plot) indeed matches the preselected one, demonstrating that 
the method in principle works. 

Fitting observed profiles. The same procedure is now used for 
fitting the observed intensity profiles (see Sect. 12. It . Although 
the MET is also clearly pronounced in Fig. |2l the best fit (not 
shown) is found just at the upper limit of the considered param- 
eter range for cr. This parameter combination is clearly unrea- 
sonable as it implies a Strehl ratio of less than 0.5 and also a 
FWHM of the PSF core that would be 30 % broader than for the 
diffraction-limited PSF. It has to be attributed to imperfections of 



the fits due to anisotropic stray-light contributions (see Sect. [Sj. 
Obviously additional constraints are required. The spatial power 
spectral density of granulation images observed with SOT would 
in principle allow to estimate the effective spatial resolution of 
the instrument. It would correspond to a line of constant FWHM 
or broadening in the y-cr-plane (e.g., a dashed line in Fig. |7ll, 
which would restrict the relevant MET and the correspond- 
ing parameter combinations to the part below. Preliminary tests 
show, however, that it is difficult to determine a strict limit for the 
spatial resolution. Here the Strehl ratio is u sed as additional con- 
strai nt. Based on pre-flig ht measurements, iTsuneta et al.l (|2008|) 
and Suemat su et al.l (|2008) state that the individual Strehl ra- 
tios of the optical telescope assembly (OTA) and the focal-plane 
package (FPP) both are 0.9 or better at a wavelength of 500 nm 
and that the comb ined Strehl ratio is thus certainly 0.8 or better 
(see Tar bell et al.L in prep., for details on FPP). Nevertheless it 
is necessary to test the optical performance under on-orbit con- 
ditions. Indeed the MET below the "knee" follows a contour of 
a Strehl ratio close to 0.8 for the Mercury transit at 555 nm (see 
Fig. IT) and 668 nm, whereas it is close to 0.85 for 450nm (see 
linear regression section of Table |4] The METes for the ecUpse 
are very similar and favour Strehl ratios of 0.84 to 0.86. 

In the following the MET above the "knee" is neglected be- 
cause of the too small Strehl ratios. Unfortunately the best pa- 
rameter set [y, cr] is less obvious than in the reference case. 
Nevertheless Fig. |7}; points at y - 0"006 and cr = 0'.'006 as 
the parameter pair that produces the best fit for that particular 
Mercury transit image. For the other cases the best pair is found 
in different locations along the MET below the "knee" and also 
the exact position of the MET varies slightly. Obviously it is dif- 
ficult to pin down a definite PSF that can reproduce all observa- 
tions. On the other hand the resulting parameter range represents 
PSFs that are very similar in terms of FWHM and Strehl ratio. 
The arithmetic averages and standard deviations of y and cr are 
listed in Table |4] for all three continuum bands for the Mercury 
transit and the total eclipse. The variation in cr can be quite sig- 
nificant, e.g. in the case of the total eclipse at 555 nm. The vari- 
ations of y and cr are not independent and should rather be re- 
placed by the deviation from the mean MET. As the MET be- 
low the "knee" follows roughly lines of constant Strehl ratio and 
broadening, the corresponding PSFs are very similar. This part 
of the MET can be approximated with a linear regression that 
accounts for all fits of the individual observations (see Table |4] 
The average and standard deviation of the slope Acr/Ay and the 
offset crc provide a range of METs, which can be translated into 
a meaningful error margin in y and cr. It is 6y ~ O'.'OOl and 
6cr ~ 0'.'002 - 0'.'003, which is of the order of the grid incre- 
ment. Owing to the larger variation in the observational intensity 
profiles for the eclipse, the METs for the eclipse at 555 nm vary 
by 6y ~ 0'.'002 and 6cr ~ O'.'OOS. 

The optimum PSF for each case can now be determined from 
the £- weighted mean value projected onto the mean MET. The 
resulting [y, cr] is essentially identical with the arithmetic aver- 
ages or at most deviates in y by just of the order of the grid in- 
crement. For the Mercury transit a cTopt = 0"008 is found for all 
three wavelengths, whereas yopt increases with wavelength from 
0"004 to 0'.'006. The results for the green and red continuum 
are very similar in Strehl ratio S (0.78-0.79) and broadening b 
of the central PSF lobe (1.4 % - 1.3 %). In contrast the blue con- 
tinuum indicates a slightly better instrument performance with 
S = 0.84 and b = 0.8 %. The intensity offset for the optimum 
PSFs is derived from the average of the individual y-cr-planes 
over all individual observations for each case at the parameter 
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Table 4. Parameters^ of the convolved non-ideal PSFs that produce the best fit of the intensity profiles across the Mercury disc and 
the total eclipse. 



450.45 nm 



555.00 nm 



668.40 nm 



450.45 nm 



555.00 nm 



668.40 nm 



Mercury transit 



Eclipse 



Convolution with Voigt function: best fit '' 


7opt 

0"opt 

b 
S 


4 (4 ± 2) 

8 (8 ± 5) 

0.3 (0.5 + 0.4) 

0.8 (0.8 ± 0.2) 

0.84 (0.84 ± 0.03) 


5 (6 ± 2) 

8 (8 ± 4) 

0.4 (0.3 ± 0.3) 

1.4(1.4 + 0.1) 

0.78 (0.78 ±0.01) 


6(6±2) 

8 (8 ±5) 

0.4 (0.4 ± 0.3) 

1.3 (1.3 ±0.1) 

0.78 (0.78 ±0.01) 


5 (5 ± 2) 

7 (7 ± 5) 

0.1 (0.2 ±0.2) 

0.9 (0.8 ± 0.2) 

0.83 (0.84 ± 0.02) 


5(4±2) 

7 (7 ± 3) 

2.5 (2.7 ± 0.4) 

0.9 (0.8 ± 0.3) 

0.84 (0.86 ± 0.04) 


4 (4 ± 1) 

7 (7 ± 3) 

3.0 (3.2 ± 0.3) 

0.8 (0.7 ± 0.2) 

0.85 (0.87 ± 0.03) 



Linear regression of MET below the "knee" 



Act /Ay 


-2.672 ±0.218 


-2.228 ±0.164 


-2.277 ±0.148 


-2.587 ± 0.309 


-2.528 ± 0.237 


-2.619 ± 0.097 


o-f 


18.2 ±1.6 


20.4 ±1.7 


21.0±1.5 


18.7 ±0.9 


17.7 ±3.2 


17.8 ±2.1 


b 


0.8 ±0.2 


1.4 ±0.3 


1.2 ±0.2 


0.9 ±0.2 


0.9 ± 0.3 


0.8 ± 0.2 


S 


0.85 ± 0.03 


0.78 ± 0.03 

r^ .. 1.,* 


0.79 ± 0.03 


0.84 ± 0.03 


0.86 ± 0.04 


0.86 ± 0.03 



7 


3±1 


4±1 


4± 1 


3±1 


3± 1 


3±1 


cr 


13±2 


12±1 


12 ± 1 


14 ±2 


12 ±2 


13 ±5 


lo 


0.4 ±0.3 


0.4 ± 0.6 


0.4 ±0.3 


0.1 ±0.2 


2.5 ± 0.4 


3.0 ±0.4 


b 


0.9 ± 0.2 


1.5 ±0.2 


1.4 ±0.2 


1.0 ±0.1 


1.1 ±0.5 


0.9 ± 0.2 


S 


0.84 ± 0.02 


0.78 ±0.01 


0.78 ± 0.02 


0.83 ± 0.02 


0.84 ± 0.05 


0.85 ± 0.04 



r 


7±1 


9±1 


9± 1 


8±1 


7±1 


7±1 


/o 


0.8 ±0.7 


0.6 ±0.3 


0.9 ±0.5 


0.1 ±0.2 


2.7 ± 0.4 


3.2 ±0.4 


b 


0.7 ±0.1 


1.3 ±0.1 


1.2±0.1 


0.8 ±0.1 


0.8 ± 0.3 


0.7 ± 0.2 


S 


0.86 ± 0.02 


0.79 ± 0.02 


0.79 ± 0.02 


0.84 ± 0.02 


0.86 ± 0.04 


0.87 ± 0.03 



7 and a in units of 10"^^ ", intensity offset /„ in %, broadening factor b of the FWHM with respect to the ideal PSF in %, and the Strehl ratio 5 . 

The values in parentheses are the arithmetic averages over the best fits for the individual observational profiles for each case. Please note that 

the errors in y and a are not independent (refer to the text for more details). The value in front of the parentheses is the £- weighted mean value 

projected onto the mean MET. 

The MET below the "knee" is approximated by (Tmedi = (Tc + (Acr/Ay) • y. The positions of the individual METs have a rms variation of 

Sy ~ O'.'OOl and 6(t ~ 0'.'002 - 0'.'003, which is of the order of the grid increment. Only the METs for the eclipse at 555 nm show a larger 

variation with Sy ~ O'.'OOl and Scr ~ 0'.'005. 

The range of reasonable PSFs is limited by the position just at the "knee" of the mean MET and on the other side by a convolution with a pure 

Lorentzian (cr = 0). 



position [yopt, cTopt]- It is only 0.3 % to 0.4 % of the mean solar 
disc intensity. 

The eclipse case indicates an excellent instrument perfor- 
mance very close to the blue continuum Mercury transit case. 
On the other hand larger intensity offsets of 2.5 % and 3.0 % are 
necessary for the green and blue continuum, respectively, clearly 
demonstrating the larger uncertainties arising from the variation 
among the eclipse profiles. In addition there could be systematic 
differences between the Mercury transit and the eclipse obser- 
vations due to the different observation dates. Some instrumen- 
tal properties change with time under space condit ions, causing 
for i nstance a gradual shift of the focus position jlchimoto et alj 
12008:) . 

Finally the range of reasonable PSFs can be limited by the 
cases just at the "knee" of the mean MET and on the other side 
by a convolution with a pure Lorentzian (cr = 0, see Sect. 14.31 ). 
The indicated instrument performance is very close to those of 
the optimum PSFs. See Table|4]for details. 

4.3. PSFs including convolution with a Lorentzian 

For the case cr -Q the deviation from the diffraction-limit essen- 
tially reduces to a pure Lorentzian with a FWHM of y. Fitting 
the observed profiles requires values for -y between 0'.'007 and 
0'.'009 (see Table Hand Fig.Hi, cf. Fig.[HJi). For the green and 
red continuum Mercury transit case the best fits imply PSFs with 



a Strehl ratio of S - 0.79 and a broadening of the central lobe 
with respect to the diffraction-limited PSF by 1.3 % and 1.2 %, 
respectively. The other cases favour Strehl ratios of 0.84 to 0.87 
and a broadening of 0.7 % to 0.8 %. The intensity offset !„ is in 
most cases somewhat higher than for the optimum PSFs for a 
convolution with a Voigt function (see Sect. 14.21 ) but stays of the 
order of 1 % for the Mercury transit and up to 3.2 + 0.4 for the 
red continuum eclipse. 



4.4. PSFs including convolution with a Gaussian 

Convolution of a diffraction-limited PSF with only a pure 
Gaussian does produce fits to the observed intensity profiles 
that are usually worse than for the analogeous cases with a 
Lorentzian (Sect. 14.31 ) or better a Voigt function (Sect. 14.21 ). The 
goodness s displayed in Fig.|8}; does not indicate any preferred 
parameter cr. Obviously a pure Gaussian is generally not a suit- 
able approximation for the non-ideal contributions to a PSF, 
especially for its wings. Although both the convolution with 
a Gaussian and Lorentzian lead to an increase in FWHM and 
Strehl ratio with their primary parameter (cr and y, respectively), 
the details of the resulting PSF profiles differ. Already for small 
values of cr the Gaussian causes the central lobe of the combined 
PSF to get so broad that the first minimum is removed and the 
central and secondary lobes essentially fuse. The higher order 
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Fig. 8. Goodness e (left column, solid lines) of the fits between 
artificial and observed intensity profiles for the Mercury transit 
2it A- 555 nm and applied intensity offset /q (right column, solid 
lines) for PSFs with convolution with a Lorentzian (top row) and 
a Gaussian (bottom row). The goodness is averaged over the in- 
dividual observations. The dashed lines are the Strehl ratio in 
the left column and the FWHM of the combined PSF in the right 
column, respectively. 



lobes are smoothed out equally fast. The Lorentzian obviously 
has the potential to produce profiles that are closer to reality. 

4.5. Application of non-ideal linearly combined PSFs 

Despite conceptual concerns the individual observed profiles 
can be fitted reasonably well when using a linear combination 
of an ideal PSF !Pdi and an additional Lorentzian (see Eq. |7] 
and Fig. |4|i. Owing to the large fluctuation among the observed 
intensity profiles (see Fig. [TJ, the best combination of a and 
j6, however, varies strongly. For instance, the observed eclipse 
profiles at 555 nm are best matched with the artificial eclipse 
experiments when adding a Lorentzian with a = 0.17 + 0.10 
and /3 = 0.07 ± 0.03, while a residuum of A/^s = (2.9 ± 1.3) % 
remains. But the range of typical parameters is extended: 
a = 0.05 - 0.30 and/3 = 0.05 - 0.10 with AI,^, = (0.5 - 4.2) %. 
Similar values are found for the eclipse at 555 nm and 668 nm 
but also for the Mercury transit at 450 nm. There are individ- 
ual cases for which a can be up to 0.35 but an increase in a is 
compensated by a relatively small value in /? and vice versa. In 
analogy to the convolution approach in Sect. 14.21 the goodness 
of fit shows this clear relation between a and /3 (see Fig. |9l). 
Obviously there is a certain ambiguity in the parameter pair 
[a,/3]. For the other Mercury transit cases at 555 nm and 668 nm 
rather values around a = 0.05 -0.10 and /? = 0.15 -0.25 
with /res - (0.2 - 1.0)% prevail. Considering all profiles for 
all wavelengths and Mercury transit and eclipse together 
the following value ranges are found: a = 0.15 (0.05 - 0.30), 
/? = 0.10 (0.05 - 0.25), (/,-es = 1.8 % (0.0 - 4.5) %. The most fre- 
quent combination of [a, yS] is [0.25,0.05] with 17 % of all cases, 
followed by [0.20, 0.05] and [0.05, 0.25] (both 13 %). 

As this approach suffers from an ambiguity between a and 
y6, additional constraints are necessary -just like the convolution 
approach in Sect. 14.21 On the other hand some cases have to 
be discarded because they are unphysical, e.g. when the Strehl 



ratio would be greater than 1 . In general /? should be equal to or 
larger than the FWHM of the diffraction-limited PSF. Hence the 
parameter domain left of the thick solid line in Fig. |9^, which 
marks /3 = FWHMji, is invalid. 



5. Discussion 

Suitable observations. Observations of the solar Umb have 
fundamental problems with respect to the determination of a 
PSF. The centre-to-limb variation induces an intensity gradi- 
ent across the image, which makes a limb observation an inho- 
mogenously illuminated and thus non-ideal test case. The corre- 
sponding stray-light level might be strongly underestimated as 
the overall light level is much smaller than at disc-centre (see 
Fig. |2]i. In addition the centre-to-limb variation complicates the 
definition of a reference intensity for the stray-light level. 

In general it is preferable to have a homogeneous back- 
ground with a well-defined proxy for the PSF. Ideal is a case 
which is not too different from regular observations. In this re- 
spect a Mercury transit seems to be the best available test case. 
The overall light level in the telescope is not significantly af- 
fected if the the FOV is much larger than the occulting Mercury 
disc. The radial geometry is very convenient as radial averaging 
offers a way to deal with the anisotropy of the PSF and to de- 
scribe it with an average intensity offset instead. A clear draw- 
back of course is that a Mercury transit does only provide the 
first 5 "of the PSF profile. 

In contrast even the far wings can be analysed for an eclipse 
image. On the other hand eclipses have a clear disadvantage if 
there are anisotropic contributions that vary across the FOV. In 
that case, as found here for SOT, the exact position of a eclipse 
terminator in the image introduces another degree of freedom. 
Also the larger relative occulted image area, which varies from 
image to image, make eclipse image depend much stronger on 
the overall light level in the telescope. 

Ideally one should try to combine the analysis of both. 
Mercury transits and eclipses, whenever possible, where images 
away from the solar limb are preferable. 

PSF models. It is not obvious which mathematical model 
is best f or representing th e non-ideal contributions of a PSF. 
Akeady iNordlundl (1 19841) made clear that PSFs of compa- 
rable quality can be constructed from a combination of two 
Lorentzians instead of Gaussians. A nother approach is used by 
iDeForest & Wills-Davevl d submitted!) , who model the PSF of the 
Transition Region and Coronal Explorer (TRACE) with the sum 
of a narrow Gaussian core and Lorentzian wings. 

An important finding is that it is indeed crucial to use real- 
istic detailed PSFs. Ideally one would assemble a reahstic PSF 
by following the optical path step by step and adding a compo- 
nent for each optical element. Unfortunately this is not possible 
in most cases as it requires detailed measurements and testing of 
each element under conditions that are very close to those of the 
observations. In a real optical system chromatic aberration re- 
sults in a systematic change of the optical properties with wave- 
length. In case of SOT, the chromatic aberration is small for the 
OTA dSuematsu et alj|2008l) and also for the BFI and can thus 
usually be neglected. An exception may be simultaneous obser- 
vations that span a large wavelengt h range, resulting in a still 
small but possibly significant effect (Ichim oto et al.l l2008'). Even 
an interchange of intensity between different wavelengths is in 
principle possible. A PSF for a certain wavelength thus could 
gain or lose energy from other wavelengths. Here the PSFs are 
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treated as purely monochromatic because the aforementioned ef- 
fects are considered negligible in comparison to other uncertain- 
ties. 

Simplifications are actually unavoidable but one should at 
least account for instrumental details such as spiders and a cen- 
tral obstruction, as it is the case for SOT. In particular the spider 
legs do not only block the telescope pupil but also scatter the 
incoming light at small angles that will inevitably broaden the 
PSF. The strong variation in eclipse profiles and the number of 
potential causes, however, limit the number of meaningful free 
model parameters. As long as an approximation produces ac- 
ceptable fits to the intensity profiles, it should thus be as simple 
as possible. For instance, the net effect of a possible defocus is in 
principle included in the PSFs determined here although it can- 
not explicitly be separated from other effects. The introduction 
of a defocus term would certainly complicate the model without 
producing better fits of the intensity profiles. 

A linear combination of a diffraction-limited PSF and a 
Lorentzian should only be used for sufficiently large j6, which 
is mostly the case when the deviations are dom inated by the see- 
ing induced by the Earth's atmosphere (see, e. g. lLangangen et aU 
l2007t iLeenaarts & Wedemever-Bohml |2005|) . For space-borne 
instruments smaller values of p and thus narrower Lorentzians 
can occur, as seen here for example of SOT. For very small y6 
the FWHM of the combined PSF gets smaller than the one for 
the diffration-limited PSF, whereas the Strehl ratio can exceed 1 . 
In those cases a linear combination is obviously unphysical. A 
Gaussian as PSF, although often used for the sake of simplicity 
and/or the lack of constraints on the PSF, is certainly not a good 
approximation, as it grossly underestimates the wing contribu- 
tions. Instead one should whenever possible prefer a convolution 
of a detailed ideal PSF and a Voigt function or possibly a more 
realistic description of the non-ideal contributions. 



cause the growth in area is more than compensated with the 
steep decrease of the PSF profile. The use of a (truncated) Voigt 
function is still justified considering that in view of the many 
sources of uncertainty i t is only a n approximation after all. 
iDeForest & Wills-Davevljsubmittedl) truncate the non-ideal con- 
tributions to the PSF wings for TRACE by applying a Gaussian 
envelope. Nevertheless a proper normalisation of the PSF re- 
quires a sufficiently large array (here a few 10"in each direc- 
tion). As a consequence the exact Strehl ratio depends on the 
extent of the PSF. Nevertheless the Strehl ratios of the PSFs con- 
structed here are of t he order of 0.8 or better and thus agree with 
iTsunetaet all (120081) . 



Stray-light estimate for SOT. Neither can one give a single 
number for the stray-light contribution but at best a rough es- 
timate for the upper limit of the stray-light contribution. As such 
one can use the difference between the observed residual inten- 
sity and the contribution due to the wings of the ideal (theoreti- 
cal) PSF. The artificial Mercury transits for the best-fitting PSFs 
(see Table Hi have intensities of 6.5 %, 10.0%, and 10.8% at 
a distance of 5" from the terminator for green, blue, and red 
continuum, respectively. These values are close to the estimated 
residual intensity in "regular" granulation images (Table|2]l. The 
small Mercury disc affects the overall light level only slightly, 
making it a good test case. Subtracting the intensities at 5" de- 
rived with the diffraction-limited PSF (see Table[3]l yields 4.7 %, 
7.8 %, and 8.2 %, which can be considered as approximate upper 
limit for the stray-light contribution at that distance for high light 
levels. The value for the green continuum is significantly higher 
than judging from the eclipse image in Fig. |3]but the latter of 
course has a lower overall light level, which causes the residual 
intensities to be smaller (see Fig.|2]i. 



Limitations of tine detailed PSFs for SOT. Anisotropic contri- 
butions and the complex dependence on the overall light level in 
the optical system and on instrumental details such as focus posi- 
tion render it simply impossible to determine only one PSF that 
is suitable for all situations. It is certainly advisable to instead 
use a range of detailed PSFs, which - as discussed in Sect. 14.21 
for the Voigt-convolution approach - is limited by the extreme 
PSFs at the MET knee and for a pure Lorentzian (see Table |4|. 
Both are still similar in terms of Strehl ratio and FWHM of the 
central lobe. Applying these PSFs in addition to the best-fit PSFs 
provide an error margin for the quantity under investigation. 

The PSFs constructed here combine the instrumental effects 
of the optical telescope assembly (OTA) and the Broadband 
Filter Imager (BFI). It is not clear to what extent the results 
can be applied to observations with the Narrowband Filtergram 
Imager (NFI) and the spectro-polarimeter (SP) because the last 
part of the optical path inside the focal plane package differs. In 
principle the PSF determination should be repeated for the other 
instruments but no suitable observations of the Mercury transit 
were found. 

A remaining uncertainty concerns the normalisation of the 
PSF to its integral, to which the far wings in principle con- 
tribute. The surface integrals of the Voigt functions diverge in 
the limit of infinite angles as the Lorentzian contribution con- 
verges too slowly. Therefore a Lorentzian needs to be trun- 
cated in order to keep the total energy of the PSF finite. Here 
it is limited by the array size of the combined PSF. From a 
practical point of view one can argue that the integrals of the 
PSFs in good approximation approach an asymptotic value be- 



Fitting procedure. Fitting a profile across a terminator is diffi- 
cult as it both has an extremely steep slope and a long tail at low 
intensities. Depending on the chosen measure either one or the 
other is weighted more. As it is not obvious which approach is 
the best, the least squares method chosen here is kept as simple 
as possible. That the choice is reasonable is supported by the fact 
that the resulting MET in the y-cr-plane indeed coincid es with 
the pre-flight measurements of the Strehl ratio (Ts unetaet al.l 
2008). Removing the denominator in Eq. (|9]l gives more weight 
to the difference between observed and synthetic profile directly 
at the steep slope. The result is a small shift of the MET in the 
y-cr-plane. For the Mercury transit the MET is shifted upwards 
by Act = O'.'OOl - 0'.'002 and to the right by Ay < O'.'OOl, which 
is just of the order of the resolution of the y-cr-plane. The Strehl 
ratio along the MET below the "knee" increases by just 0.01 
to 0.03. The quality measure given in Eq. (|9]l produces balanced 
fits, whereas removing the denominator can cause significant de- 
viations in the tail. The eclipse is slightly more susceptible to 
changes of the quality measure e, producing shifts of the MET 
of Ay = O'.'OOl - 0'.'004 and Act = 0'.'002 - 0'.'008 and a change 
of the Strehl ratio of up to ~ 0.10 in some cases. Reducing the 
spatial range that is considered for e from 10" to 5" effects the 
position of the MET much less with shifts comparable to the 
above-mentioned uncertainties for the Mercury transit. 

In summary it can be concluded that the Mercury transit case 
is less susceptible to details of the fitting procedure than the 
eclipse case. 
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Intensity offset. Although the non-ideal PSFs provide a much 
better fit of the observed intensity profiles than the ideal ones 
small quantitative discrepancies remain. The remaining differ- 
ences reflect the fact that the detailed PSF models constructed 
here are still a necessary simplification whereas the real PSFs 
additionally include anisotropy of the residual light and thus de- 
pend on the position in the FOV. The latter would require more 
constraints, which are usually - as in the present case - not avail- 
able. The introduction of the constant intensity offset Iq offers a 
way of capturing the net effect of the position-dependence and 
improves the fit of the intensity profiles. Furthermore !„ absorbs 
uncertainties in the determination of reference intensities and of 
the fitting procedure but also potential small contributions from 
the very far PSF wings. The values for /q, however, are mostly 
only of the order of 1 % of the mean solar disc intensity or less 
for the Mercury transit (see, e.g., Fig.|7};). The variation of the 
mean count rate and the anisotropy of the stray-light produce 
a larger remaining uncertainty, requiring intensity offsets of up 
to ~ 3.6 % (for the optimum PSFs with convolution of a Voigt 
function). 



Anisotropic contributions. For an observed profile anisotropic 
contributions due to instrumental scattered light - although small 
- are expected and indeed found in case of SOT (Fig. [3]). Adding 
an intensity contribution of the form Iadd(x) - x j-J to the refer- 
ence case changes the slope of the synthetic intensity profiles 
across the terminator and removes the otherwise clear fading 
of the MET on the upper branch. The parameter pair produc- 
ing the best fit, i.e. the minimum error, in the j-cr plane moves 
along the MET to the upper right to Strehl ratios that are cer- 
tainly unreasonable. In reality the slope changes are certainly 
more complex than the simple linear contribution tested here. 
Consequently the used synthetic intensity profiles cannot per- 
fectly match the observed profiles. There will always be a re- 
maining (small) discrepancy. An exact match would require a 
detailed model of the anisotropic contributions, which would de- 
mand for more observational constraints on the PSF properties. 
Obviously anisotropic contributions severly complicate the de- 
termination of a reasonable PSF. 



Alternative methods. The direct determination of a (two- 
dimensional) PSF from an observational image of an eclipse 
would in principle be superior to the forward approach presented 
here. In most cases, however, a direct determination is rather 
hopeless in view of uncertainties that arise from the PSF wings 
but also from detector sensitivity, data sampling and many more. 
All these effects make it difficult to retrieve a definite analytical 
fit of a measured intensity profile. While a radially symmetric 
PSF can be derived this way but asymmetric contributions due 
to, e.g., instrumental stray-light are usually hard to determine. 
Obviously non-radial features such as the imprint of spider legs 
cannot be captured this way. 

One has to keep in mind that the situation is quite different 
from night time observations for which stars, i.e. point sources, 
offer an easy way to determine a PSF. It may be possible to di- 
rectly measure a two-dimensional PSF for a solar telescope by 
taking images of background stars far away from the solar disc 
but there the conditions are fundamentally different from those 
of regular solar observations. As the dependence on the overall 
light level in the telescope effects in particular the PSF wings, 
a measurement outside the solar disc may produce a PSF that 
cannot be applied to observations of the bright solar disc. A po- 
tential technical problem is that solar telescopes and in particular 



the detectors are usually not designed for such low-intensity sit- 
uations so that the low signal-to-noise ratio prevents the direct 
measurement of a PSF. 

The analysis of very small features in principle provides 
an indicator of the effective spatial resolution and with it the 
FWHM of the PSF. ISuematsu elal] ([2008) use a intensity pro- 
file of a small G-band bright point and conclude that the OTA 
of S OT indeed oper ates close to the diffraction-limit (see also 
Ichi moto et al.ll20()8l) . Unfortunately this method does not yield 
further information about the detailed structure of the PSF. The 
same is true for the spatial power spectral density of quiet Sun 
granulation images. 

Another common method is to determine a PSF by matching 
the intensity distribution based on numerical simulations of the 
solar atmosphere with observations. This approach should only 
be considered in the total absence of other independent methods 
as it assumes the numerical simulations to be perfectly realistic. 

6. Conclusion 

The analysis of Mercury transit and eclipse broadband filter- 
grams shows that SOT onboard the Hinode satellite is an ex- 
cellent instrument, which performs close to the diffraction limit. 
In this work PSFs are constructed by a convolution of the ideal 
diffraction-limited PSF with a narrow Voigt function, which ac- 
counts for non-ideal contributions due to instrumental stray- 
light. The comparison of synthetic and observed intensity pro- 
files across the Mercury and lunar terminator does unfortunately 
not produce a unique solution for the best-fitting PSF, which can 
account for all observational conditions. Likely reasons are the 
anisotropy of the residual light and the dependence on the over- 
all light level. A realistic PSF would even depend on the position 
on the solar disc, owing to the change of the overall light level 
in the telescope due to the centre-to-limb variation of the contin- 
uum intensity. Obviously it is impossible to provide a stray-light 
estimate as just a single (wavelength-dependent) number. 

Instead of a single PSF for each wavelength one should 
rather use a range of PSFs, which allow to translate the uncer- 
tainty in PSF properties into an error margin for the quantity un- 
der investigation. Fortunately the differences between synthetic 
and observed intensity profiles form a distinct relation of the 
Voigt parameters y and cr, which is named here the "minimum 
error trench" (MET). The empirically estimated Strehl ratio of 
the optical system serves as additional constraint, limiting the 
range of reasonable PSFs to a part of the MET where the PSFs 
are very similar in terms of Strehl ratio, FWHM of the central 
lobe, and thus their effect on the intensity distribution. There 
are remaining uncertainties due to the anisotropy and light-level- 
dependence of the residual intensity, which result in an intensity 
offset during the profile fitting process. For the Mercury transit 
the offsets are usually only of the order of 1 % of the solar disc 
intensity or less, whereas the eclipse images cause values of up 
3.6 % with a pronounced variation among the individual images. 
Obviously the Mercury transit is the more reliable test case. 

Nevertheless reasonable PSF estimates can be constructed 
(see Table nil. They all require only narrow non-ideal contribu- 
tions, indicating an excellent performance of the instrument. 
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Fig. 1. Hinode/SOT observations from top to bottom: a) limb observation, c) the total eclipse from March 19th, 2007, e) the partial 
eclipse from February 17th, 2007, and g) the Mercury transit from November 8th 2006 on the solar disc (x - -52 ", y = -426 ", 
ju - 0.9) and i) at the solar limb. For all cases an exemplary image for a wavelength of 555 nm is selected, except for the partial 
eclipse which was only observed in G-band (431 nm). The panels to the right show intensity profiles of the image shown on the 
left (black solid) and for all other available images of the same kind (dark-grey). The profiles are averaged over all rows, aligned 
perpendicular to the terminator/limb. The intensity is normalised to the intensity of the bright solar disc close to the limb / terminator. 
The left half of the panels h and j, show averages over the three rows and columns through the centre of the Mercury disc, whereas 
the radial averages can be seen in the right half. The dotted lines are for orientation only. For comparison profiles from the total 
eclipse observed in G-band are plotted in panel f (dashed lines). 
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Fig. 2. Residual intensity as function of mean 
count rate at a distance of 5 " from the (eclipse) 
terminator, which corresponds to the centre 
of the Mercury disc. The panels in the up- 
per row (a-c and the leftmost column (d, g, j) 
show the residual intensity /^s,^ as function of 
the mean count rate C for the individual im- 
ages of the Mercury transit (circles), the to- 
tal eclipse (squares), the partial eclipse (dia- 
monds, 431 nm only), and the solar limb (trian- 
gles) at all six BFI bands. Symbols are filled for 
yu > 0.5 and unfilled and dark-grey otherwise. 
The adopted /i-value for a data point is the aver- 
age over the whole individual image. A image 
of the solar limb thus can have fi < 0.0. A rela- 
tion between residual intensity and mean count 
rate is derived by linear regression (solid lines), 
where the solar limb data is neglected. It is jus- 
tified by the uncertainties in choosing the right 
reference intensity (see text for details). More 
information is shown for the three wavelengths 
450 nm, 555 nm, and 668 nm. First the helio- 
centric position fi is plotted as function of C in 
the middle column (e, h, k but not b) for the 
eclipse and Mercury transit data and for images 
of "regular" quiet Sun (dots). The relation be- 
tween /res.i and C is then used to estimate I^^s.a 
from the mean count rate of the regular (unoc- 
culted) images. The results are presented in the 
right column (f, 1, 1). As abscissa fi is chosen. 
Finally a linear regression is performed for the 
data points of regular granulation (solid lines). 
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Fig. 6. Goodness s of the fits for the reference case of the Mercury transit at A - 555 nm. It represents the difference between a 
prescribed synthetic intensity profile for a parameter pair of y = 0'.'005 and cr - O'.'Ol and all other profiles. They are derived with 
non-ideal PSFs calculated by convolution with a Voigt function, a) The goodness s as defined in Eq. (|9]l as function of j and cr, b) 
close-up of the lower left corner, c) s along the "minimum error trench" (MET, marked with white dotted line in panel a), and d) 
the applied intensity offset I^ along the MET. Overplotted in panel a and b are contours of constant Strehl ratio (solid, black and 
white) and of constant broadening (dashed) of the combined PSF. The broadening is given in percent points increase of the FWHM 
of the PSF with respect to the diffraction-limited PSF. The circles in panel a and b mark the parameter pair with the smallest e below 
the "knee" of the MET, while the "knee" is indicated by a triangle. The dot-dashed lines in panel a are for constant Voigt damping 
parameter a, whereas the dotted vertical lines in panels c and d show the rough location of the MET "knee". 
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Fig. 7. Goodness s of the fits between artificial and observed intensity profiles for a selected image of the Mercury transit at 
A = 555 nm. The synthetic profiles are generated with non-ideal PSF that are calculated by convolution with a Voigt function. See 
Fig.|6]for a detailed explanation of the figure. 
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Sven Wedemeyer-Bohm: Point spread functions for SOT 
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Fig. 9. Goodness e of the fits between artificial 
and observed intensity profiles for the Mercury 
transit al A = 555 nm (cf. Fig. [6}. The non- 
ideal PSF contribution are calculated by lin- 
ear combination with a a) The goodness e as 
defined in Eq. ^ as function of a and /3, b) 
E along the "minimum error trench" (MET, 
marked with white dotted line in panel a), and 
c) the applied intensity offset !„ along the MET. 
Overplotted in panel a are contours of constant 
Strehl ratio (solid) and of constant broaden- 
ing (dashed) of the combined PSF. The broad- 
ening is given in percent points increase of 
the FWHM of the PSF (FWHMpsf) with re- 
spect to the one of the diffraction-limited PSF 
(FWHMji). All parameter combinations with a 
P smaller than FWHMji (with except for com- 
binations with a = 0) produce unphysical PSFs 
with FWHMpsF < FWRMj,. The cases are 
found left of the thick solid line, which marks 
FWHMpsF = FWHMj, (in all panels). 



